Table 1 fcont') 

/* 

* strip path or prefix from pn, return len: pr_alignO 
*/ 

static 

5 stripname{pn) stripname 

char *pn; /* file name (may be path) */ 



{ 



register char *px, *py; 



10 py = 0; 

for (px = pn; *px; px++) 
if(*p X == •/•) 

py = px + 1; 

if(P7) 

15 (void) strcpy(pn, py); 

return(strlen(pn)); 

.= > 
: 20 



35 
40 
45 
50 
55 
60 
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Table 1 (contM 

* cleanupO - cleanup any trap file 

* getseqO — read in seq, set dna, len, maxlen 

* g_calloc0 ~ callocO with error checkin 

5 * readjmpsO ~ get the good jmps, from tmp file if necessary 

* writejmpsO ~ write a filled array of jmps to a tmp file: nwO 

include "nw.h" 
#include <sys/file.h> 

10 

char *jname = Vtmp/homgXXXXXX"; /* tmp file for imps */ 

FILE *fj; 

int cleanupO; /* cleanup tmp file */ 

15 long lseekO; 

/* 

* remove any tmp file if we blow 
*/ 

20 cleanup® cleanup 

int i; 

{ 

if(S) 

(void) unlinkQname); 

25 exit(i); 
i 

I* 

* read, return ptr to seq, set dna, len, maxlen 
30 * skip lines starting with ' < or ' > * 

* seq in upper or lower case 
*/ 

getseq(file, len) getseq 
-*— /*filen 



{ 



*len; /* seq len */ 



line[1024I, *pseq; 
*px, *py; 
int nalgc, tlen; 

FILE *fp; 

if (ft> = fopen(fiIe, , •r , ')) == 0) { 

fprintf(stderr,"%s: can't read %s\n", prog, file); 
exit(l); 

} 

tlen = natgc = 0; 

while (fgetsOine, 1024, fp)) { 

if (*line == •;' 1 1 *Une == '<• 1 1 *line == •>') 



for (px = line; *px ! = '\n'; px+ +) 

if (isupper(*px) 1 1 islower(*px)) 
tlen++; 

} 

if ((pseq = malloc(Cimsigned)(tlen+6))) == 0) { 

fprintf(stderr,"%s: mallocO failed to get %d bytes for %s\n\ prog, tIen+6, file): 
exit(l); 

} 

pseq[0] = pseq[l] = pseq[2] = pseq[3] = '\0'; 
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Table 1 (conf) 

py - pseq + 4; 
*Ien = tlen; 
rewind(fp); 

while (fgetsOine, 1024, fp)) { 

if (*Iine == ';* || *line == '<• || *line == '>*) 
continue; 

for (px = line; *px != *\n'; px++) { 
if (isupper(*px)) 

*py++ = *px; 
else if <islower{*px)> 

*py++ = toupper(*px); 
if (index("ATGCU",*(py-l))) 

natgc++; 

} 

} 

*py ++ = <\Q'; 



If J > v > 

20 (void) fclose(fp); 

dna = natgc > (tlen/3); 
return(pseq+4); 

} 

25 char * 

g_caIloc(msg, nx, sz) gcalloc 
char *msg; /* program, calling routine */ 

int nx, sz; /* number and size of elements */ 

char *px, *calloc0; 



if ((px = calloc((unsigned)nx, (unsigned)sz)) = = 0) { 
if(*msg){ 

n - lprintfiCstderr, " %s: g_calloc0 foiled %s (n= %d, sz= %d)\n", prog, msg, nx, sz); 

» exit(l); 



} 
} 

return(px); 



* get final jmps from dxQ or tmp file, set ppO, reset dmax: mainO 

readjmpsO readjmps 



siz, iO, il; 



• i, j, xx; 
50 if(rj){ 

(void) fclose(fj); 

if ((fd = open(jname, O RDONLY, 0)) < 0) { 

fprintf(stderr, " %s: can't openO %s\n", prog, jname); 
cleanup(l); 

55 } 
} 

for (i = iO = il = 0, dmaxO = dmax, xx = lenO; ; i++) { 
while (1) { 

for (j = dx[dmax].ijmp; j > = 0 && dx[dmax].jp.xDl > = xx; j-) 
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